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ABSTRACT 

In the standard model of cosmology, the Universe is static in comoving coordinates; expan¬ 
sion occurs homogeneously and is represented by a global scale factor. The baryon acoustic 
oscillation (BAO) peak location is a statistical tracer that represents, in the standard model, a 
fixed comoving-length standard ruler. Recent gravitational collapse should modify the metric, 
rendering the effective scale factor, and thus the BAO standard ruler, spatially inhomogeneous. 
Using the Sloan Digital Sky Survey, we show to high significance (P < 0.001) that the spatial 
compression of the BAO peak location increases as the spatial paths’ overlap with superclus¬ 
ters increases. Detailed observational and theoretical calibration of this BAO peak location 
environment dependence will be needed when interpreting the next decade’s cosmological 
surveys. 
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1 INTRODUCTION 


The choice of a spacetime coordinate system for the Universe that 
enables the expansion to be represented via a global scale fac- 
tor a as a func tion of just one coordinate (cosmological time t, 
feemaitrei e.g. 1 19271) is extremely convenient. On large enough co¬ 
moving length scales, statistical spatial patterns are fixed in this co¬ 
ordinate system. Thus, the theory of primordial de nsity fl uctuations 
leads to that of baryon acoustic oscillations feisenstein & Hti 19981 : 
BAOs). The BAO peak in the two-point spatial auto-correlation 
function £ was c learly detected in the Sloan Digital Sky Survey 
feisenste in et alj| 2005f) an d the Two-Degree Field Galaxy Redshift 
Survey fcole et al l 1 2005) : earlier surveys may have detected this 
too, lEinasto et al.l 1 19970 . BAOs now constitute one of the most 
important tools for making cosmological geometry measurements, 
especially for upcoming observational projects such as the space 
mission Euclid I Refregieret alj2010l) and the ground-based instru¬ 
ments DESI feevi et a lJl2013h . 4MOST t ide Jong et alJl2012h . and 
the LSST l lTvson et al.1 l2003h . The BAO peak location of about 
105/r 1 Mpc (where h is the Hubble constant Hq expressed in units 
of 100 km/s/Mpc) is commonly expected to be a large enough co¬ 
moving length scale for it to provide a fixed comoving ruler in the 
real Universe. 


* BFR: during visiting lectureship; JJO: during long-term visit. 


However, the validity of the BAO peak location as a stan¬ 
dard ruler depends on the validity of the assumed cosmologi¬ 
cal metric (differential rule for measuring lengths) in the context 
of the real Universe, which is lumpy (cf. the “fitting problem”, 
fellis & StoegeJ 19871) . Scalar averaging is a general-relativistic for¬ 
malism that extends beyond the standard cosmological model, by 
allowing a spatial section of the Universe at a given time to have an 
inhomogeneous metric and calculatin g background- free volume- 
weighted averages of scalar variables feuchertll200ll . [20081 ). The 
univariate scale factor a(t) is replaced by an effective, environment- 
dependent volume-based scale factor ao(t) oc V^ 3 , dependent on 
both the choice of compact spatial domain D of volume Vo and on 
time. Without this extension, the cosmological, comoving metric 
is forced (by definition) to be rigid in comoving coordinates, i.e. 
inhomogeneities in the matter distr ibution are not allowed to “tell 
comoving space how to curve” (e.g. lBuchert & Carforal2008h . Ap¬ 
plying scalar averaging to an observationally standard power spec¬ 
trum that statistically represents density fluctuations at an early 
epoch implies that even for spatial domains as large as the BAO 
peak length scale, the environment dependence of the scale fac¬ 
tor should be observationally detectable, i.e. a m < fle is expected, 
where M. (“Massive”) and £ (“Empty”) represent overdense and 
underdense spatial regions, respectively. For example, adopting lcr 
initially overdense (yVt) and underdense (£) fluctuations in a spher¬ 
ical domain of diameter « 105/r 1 Mpc and using equations (2), 
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Figure 1. Compression of the baryon acoustic oscillation peak. Upper 
panel: BAO peak (cubic-subtracted correlation function) for pairs of lu¬ 
minous red galaxies whose paths either overlap with superclusters by 
co > o) m i n = 60 h~ [ Mpc or are entirely co ntained with in th e sup erclusters. 
The overlap co is the chord length defined iRoukema et alj|2015l Sect. 2.3, 
Fig. 1) by the intersection of the path joining two LRGs and the superclus¬ 
ter modelled as a sphere. Individual curves represent 32 bootstrap resam¬ 
plings of the supercluster catalogue (235 objects) and the “random” galax¬ 
ies (484,352 selected from 1,521,736). The real galaxies (30,272) are not 
resampled. Most of the curves peak sharply at 95 h~ l Mpc; a few peak at 
85/i _1 Mpc. The high amplitude (in comparison with the lower panel) is 
consistent with biasing that modifies the amplitude of £. Lower panel: BAO 
peak for the complementary subset of galaxy pairs. The peak occurs at the 
standard value of about 105/z -1 Mpc. 


(13), (32), (50), and (54) of iBuchert. Navet & Wiegaiidl J2013h to 
integrate the Raychaudhuri equation (9) of the same paper gives a 
relativistic Zel’dovich approximation estimate of 


a Ml a & ~ 0.91 . (1) 

In other words, in the scalar averaging approach, one way in which 
matter inhomogeneities are expected to affect the large-scale geom¬ 
etry and dynamics is to shrink the curved-space volume of overden¬ 
sities on the BAO scale by about 1 - (a m,/ as) 3 ~ 24%, leading to an 
expected shift in the BAO peak location to a lower scale by some¬ 
what below or above 9% for galaxy pairs that weakly or strongly, 
respectively, overlap with typical BAO-scale overdensities. 



Minimum overlap co mjn (Mpc/h) 


Figure 2. Overlap dependence of the BAO peak shift. BAO peak shift 
A^ := ^non-sc — -Ssc, where ,s’ sc and Snon-sc are the median estimates of the 
centres of the best-fit Gaussians to £ — £3 for LRG pairs that overlap su¬ 
perclusters (sc) and those that do not (non-sc), respectively. The error bars 
show a robust estimate of the standard deviation, cr(A.s), defined here as 
1.4826 times the median absolute deviation of As. At each cu m i n , these 
statistics are calculated over 32 bootstrap resamplings of the observational 
data. Four out of the 320 Gaussian fits in the supercluster-overlap case failed 
and were ignored in calculating these statistics; no failures occurred for the 
320 non-supercluster-overlap cases. A linear least-squares best fit relation 
As = 4.3 h~ [ Mpc+0.07w min is shown with a red line. The most signifi¬ 
cant individual rejection of a zero shift is As = (6.3 ± 2.0)/z —1 Mpc for 
co m j n = 30 h~ [ Mpc, a 3.1cr (Gaussian) rejection. Since bootstraps are used, 
this estimate is c onservative: o~(A^) is expected to be an overestimate of the 
true uncertainty dpishcr et al.ll 1994 Sect 2.2). A 9% shift [Eq. [Q] would 
give As = 0.09m. Since a> > m nnni a scalar-averaging lower expected limit 
As > 0.09w min is shown as a green line. 


An environment-dependent effect has recently been detected 
as a six percent compression of the BAO peak location for spatial 
paths that touch or overlap superclusters of luminous red galax¬ 
ies (LRGs) in the Sloan Digital Sky Survey dRoukema et TlboiSl: 
SDSS; environment dependencejtf £ at smaller scales has also been 
detected in the SDSS. Ichiang et all2015l) . 

In this Letter , we check whether the compression is dependent 
on the minimum overlap between spatial paths and superclusters, 
as it should be if the effect is induced by the statistically overdense 
nature of the superclusters. 


2 METHOD 


We modify the previous method dRoukema et al.ll7oi5L Sect 2) in 
order to allow stronger overlaps. As in the original method, we cal¬ 
culate the correlation function f of the “bright” sample of LRGs 
in the SDSS Data Release 7 (DR7) for pairs of LRGs selected for 
overlap IRouke ma et alJl2015l Sect. 2.3, Fig. 1) (or non-ov erlap) 
of superclusters in the survey IjNadathur & Hotchkis s! 120141) . us¬ 
ing the Landy & Szalay estimator (Land^_<&S!zalav 1993) on real 
and “random” (artificial) catalogues dKazin et al ■ 20 hT Comoving 
separations s are ca l culated assumin g the standard ACDM model 
dSpergel et alj|2003l : [Ade et alj|20l4) with matter density parame¬ 
ter D m0 = 0.32 and dark energy parameter Q A0 = 0.68. A best-fit 
cubic £ 3 (s) over separations 5 < 70 hr 1 Mpc and s > 140 hr 1 Mpc 
(i.e. excluding the peak), is found for the tangential signal (pairs 
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< 45° from the sky plane). The proce dure is repeated, bootstrap re - 
sampling the supercluster catalogue dNadathur & Hotchkissll2014l) 
and the “random” LRG catalogue, several times. The BAO peak 
location is estimated from the medians and median absolute de¬ 
viations of the centres of the best-fit Gaussians to %(s) - ^(j). 
In this work, minimum overlaps oi min in the range 10/r 1 Mpc 

< w m i n < 100/r 1 Mpc rather than aw = I/7 1 Mpc are consid¬ 
ered. In order that £ be defined for i < aw for these high values 
of ow, we consider a pair of LRGs joined by a comoving spatial 
path entirely contained within a supercluster to satisfy the overlap 
criterion. 


3 RESULTS 

Figure Q] shows that for a minimum overlap oi m j n = 60 h 1 Mpc, 
the BAO peak is shifted to lower separations v (panel a) than 
for the complementary set of LRG pairs (panel b). The shift is 
clearer than for the earlier analysis, which had aw = 1/t -1 Mpc 
dRoukema et al ][2~015| . Fig. 8). Requiring a stronger overlap yields 
a stronger shift. 

Figure [2] shows the dependence of the shift As on ow. 
The BAO peak shift for supercluster-overlapping LRG pairs ap¬ 
pears to increase from As ~ 6 —7/i 1 Mpc for aw i 50/7 1 Mpc to 
As k1 1 lr' Mpc for greater overlaps. The Pearson product-moment 
correlation coefficient of As and aw is 0.87, with a probability of 
P ~ 0.0008, i.e. a positive correlation is detected to high signif¬ 
icance. Thus, the environment dependence of the BAO peak shift 
is confirmed. Unfortunately, the uncertainties shown in Fig. [2] are 
too high to infer the details of this correlation from the present 
data and analysis. The slope and zeropoint of the linear best fit are 
0.073 ± 0.040 and 4.3 ± 2.0/? 1 Mpc, respectively. 

As a rough guide to what is expected from scalar averaging, 
we can use the 9% shift estimate from Eq. (jTji, which would give 
As = 0.09ai, where the overlap path lengths are approximated as 
corresponding to lcr overdense regions on the BAO scale, even 
though in reality, the overdense regions are superclusters, some 
smaller and some larger than this scale. Since oi > ow, this implies 
a rough scalar-averaging lower expected limit of As > 0.09oi min , 
shown as a green line in Fig. [2] 


4 DISCUSSION 


Since the BAO peak location serves as a major tool for cosmo¬ 
logical geometry measurements, it is clear that its environment 
dependence will need to be observationally calibrated and cor¬ 
rectly modelled theoretically. It is possible that the effect could 
also be interpreted within the standard ACDM model, as is the 
case for many large-scale phenomena. For ex ampl e, observed 
super voids on the 200-300/7“* Mpc scale aNadathur & Hotchkisa 
120141 : [Szagudi_etal, 20151) ca n be interpreted within the ACDM 
model IHotchkiss et al.112 015l). al t hough their occurrence is ex¬ 
pected to be rare l lSzapudi et al.112015h . In contrast, the study 
of SDSS DR7 “dim” (or "bright”) LRGs via Minkowski func¬ 
tionals on scales ranging up to the BAO peak scale, within a 
500/7“' Mpc (or 700/7“ 1 Mpc, respectively) diameter region, shows 
3-5.5cr (or 0.5-2.5<x) inconsistencies with ACDM simulations 
dWiegand, Buchert & Qstermannll2014l . Table 1). Minkowski func¬ 
tionals have more statistical power than lower order statistics that 
are commonly used in analysis of large-scale structure, such as 


the two- and three-point correlation functions, the correlation di¬ 
mension or percolation (“friends-of-friends”) analyses. This is be¬ 
cause all the 77 -point correlation functions would be needed in 
order to represent the statistical geometrical information that the 
Minkowski functionals contain. 

A possible avenue to studying the environment-dependent 
BAO shift within the standard approach would be to carry out a 
Fourier analysis rather than using the two-point correlation func¬ 
tion. This would require the development of a supercluster-overlap- 
dependent Fourier analysis method. Another alternative, to avoid 
having to determine the position of the peak itself, would be com¬ 
parison of radial to tangential correlation functions directly. This 
would require correcting for peculiar velocity effects, which are 
highly anisotropic with respect to the observer. 

Interpreting the environment dependence of the BAO peak lo¬ 
cation reported in this Letter within the standard ACDM model 
would require the comoving length scale at which the Universe is 
rigid in comoving coordinates to be pushed up to a scale greater 
than 105/7“* Mpc. The environment dependence (e.g. the Aj(oi min ) 
relation) would have to be modelled within a rigid comoving back¬ 
ground that can only exist at larger scales, at which no sharp sta¬ 
tistical feature that can function as a standard ruler is presently 
known. This leads to a Mach’s principle type of concern that at 
recent epochs, it is difficult to have confidence that the standard 
comoving coordinate system is correctly attached to an observa¬ 
tional extragalactic cat alogue (peculia r velocity flow analyses indi¬ 
cate similar concerns. [Wiltshire et al]|2013l ). Interpretation within 
the scalar averaging approach should be easier because its descrip¬ 
tion of fluctuation properties and the cosmological expansion rate 
is environment-dependent. 

Nevertheless, within the rigid comoving background frame¬ 
work (i.e. the standard model), small shifts in the BAO peak lo¬ 
cation have been predicted analytically and from A-body simula- 
tions l lDesiacaues et al.ll20ld : [sherwin^Za ldarriaga 2012b, while 


BAO reconstruction techniques 1 Padmanabhan & White 1 20091 : 
Padmanabhan e t all2012l :ISchmit tfull e t al.l20lH )~ have been devel¬ 
oped to attempt to evolve galaxies’ positions backwards in cosmo¬ 
logical time, using a blend of theoretical calculations and A-body 
models. The expected mean shift in the BAO peak location is less 
than one percent, i.e. an order of magnitude less than what we find 
for the shift conditioned on oi m i n > 6 O /7 1 Mpc. The amplitude 
of the shifts found in these calculations is constrained by the as¬ 
sumption that curvature averages out on the assumed background, 
i.e. that a conservation la w for instrinsic curvature holds globally 
dBuchert & Carforall2008h . 

Theoretical work is underway in the scalar averaging ap¬ 
proach, which general-relativistically extends the standard model, 
allowing the restrictive assumptio n of a conservation law for in¬ 
trinsic curvature to be dropped dBuchert & Carforal 2008). The 
physical origin of curvature deviations from the background on 
scales as large as the BAO scale can then be thought of as fol¬ 
lowing from the non-existence of a conservation law for intrin¬ 
sic curvature. To reconstruct the primordial comoving galaxy po¬ 
sitions more accurately than in the standard model, i.e. to allow 
flexible comoving curvature that varies with the matter density 
and the extrinsic curvature tensor acr oss a spatial slice, relativis¬ 
tic Lagrangian perturbation theory ([B uchert & Oste nnanr] 120121 : 
iBuchert. Navet & Wiegandll2013l : lAlles et all2015b is available for 
analytically guided calculations. A-body simulations in which the 
growth of inhomogeneities is matched by inhomogeneous metric 
evolution will most likely also be needed to develop numerical con¬ 
fidence in what could be called “relativistic BAO reconstruction”. 
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5 CONCLUSION 

No matter which approach is chosen, analytical, numerical and ob¬ 
servational work will be required if the BAO peak location is to 
correctly function as a standard ruler for cosmological geometrical 
measurements, since the evidence is strong (P < 0.001) that it is 
strongly affected by structure formation. Moreover, the formation 
of superclu sters— in reality, filamentary and spiderlike distributions 
of galaxies dEinasto et al.l2014l) rather than the spherically symmet¬ 
ric objects assumed here for calculational speed—can now be tied 
directly to a sharp statistical feature of the primordial pattern of 
density perturbations. 
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